Spontaneous phase oscillation induced by inertia and time delay 
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We consider a system of coupled oscillators with finite inertia and time-delayed interaction, and 
investigate the interplay between inertia and delay both analytically and numerically. The phase 
velocity of the system is examined; revealed in numerical simulations is emergence of spontaneous 
phase oscillation without external driving, which turns out to be in good agreement with analytical 
results derived in the strong-coupling limit. Such self-oscillation is found to suppress synchronization 
and its frequency is observed to decrease with inertia and delay. We obtain the phase diagram, which 
displays oscillatory and stationary phases in the appropriate regions of the parameters. 
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Systems of coupled nonlinear oscillators, serving as a 
prototype model for various oscillatory systems in nature, 
have been known to exhibit remarkable phenomena of 
collective synchronization, which have been observed in 
a variety of physical, biological, and chemical systems [Q . 
Among them there are systems where finite inertia and 
time delay are relevant. For example, in a supercon- 
ducting junction network the capacitance plays the 
role of inertia while time delay naturally arises in some 
physical as well as biological systems, where finite time 
interval is actually required for the transmission of in- 
formation 1^. This has motivated recent studies of the 
effects of finite inertia and of time delay on synchroniza- 
tion: In a system with inertia suppression of synchroniza- 
tion and emergence of hysteresis has been observed [Q, 
whereas a multitude of coherent states with different syn- 
chronization frequencies and suppression of the collective 
frequency have been reported in systems with delay 
Note, however, that the effects of the two have been con- 
sidered separately and independently of each other. The 
interplay between the two has not been explored in the 
most general system with both inertia and time delay. 
For a single oscillator with inertia, it is known that de- 
layed restoring force can destabilize a fixed point and 
produce oscillation |^ ; this raises an interesting question 
as to how such oscillatory behavior affects collective syn- 
chronization in a set of globally coupled oscillators with 
inertia and delay. 

This paper is to investigate the interplay between in- 
ertia and retarded interaction and their combined effects 
on synchronization. We consider a system of globally 
coupled oscillators, each possessing finite inertia and in- 
teracting with others via time-delayed interactions. To 
explore the interplay, we first examine the temporal be- 
havior of the phase velocity, and obtain a nonvanishing 
ac component in the absence of external periodic driv- 
ing. Such spontaneous oscillation is found to suppress 
synchronization, and its frequency observed to diminish 
with inertia and delay. 

We begin with the set of equations of motion governing 
the dynamics of N coupled oscillators, the zth of which 



is described by its phase 0i (i = 1, 2, TV): 
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where /i represents the magnitude of (rotational) inertia 
relative to damping. The third term on the left-hand side 
denotes the global coupling of strength K/N between 
oscillators, indicating that each oscillator interacts with 
others only after the retardation time r. The term uji 
on the right-hand side denotes the intrinsic frequency 
of the ith oscillator, and is randomly distributed over 
the whole oscillators according to g{Lu), which is assumed 
to be symmetric about uj — 0, and concave aA, lo — 0. 
In the absence of time delay and inertia (r = /i = 0), 
Eq. (U) exactly reduces to the Kuramoto model 0, for 
which analytical results are available. The system with 
either time delay or inertia (but not both) also exhibits 
synchronization for the coupling strength strong enough. 
The phases are stationary or monotonic in time 
indicating that the system either relaxes to the minimum- 
energy configuration or displays constant phase velocity. 

Here we investigate the dynamics of phases in the pres- 
ence of both inertia and delay, first by means of numer- 
ical simulations. For convenience, we take the Gaussian 
distribution with unit variance (a'^ = 1) for gluj) and 
choose the coupling strength K = 3, which is larger than 
the known critical values [^,1). We thus probe the in- 
terplay between inertia and delay in the synchronized 
state. Equation (|l|) has been integrated with discrete 
time steps of St — 0.01. In computing the order parame- 
ter, Nt = 10^ time steps have been used at each run, with 
the data from the first 9.6 x 10^ steps discarded. We have 
varied both 6t and Nt to verify that the steady state has 
been attained and performed twenty independent runs 
with different initial configurations, over which averages 
have been taken. In this manner, we have computed the 
power spectrum of the phase velocity: 
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where = / 4'j{t)e'^'^'^^*dt is the Fourier component 

of the phase velocity at frequency /. 

Figure 1 shows the obtained power spectrum in a 
system of iV = 100 oscillators, for coupling strength 
if = 3 and time delay r = 2 at various inertia values, 
H = 0.4, 0.6, and 0.8. (We have considered the size N up 
to 3200, and confirmed the absence of appreciable finite- 
size effects for N > 100). For small inertia (/i ^ 0.3), 
the power spectrum possesses only the zero-frequency 
component. Namely, the system displays stationary or 
monotonic behavior of the phase, like the one with de- 
lay only Q. When the inertia is raised beyond 0.3, on 
the other hand, several peaks develop at nonzero fun- 
damental and harmonic frequencies, and the oscillatory 
phase emerges. As /i is increased further, those peaks 
shift to the lower frequency sides, eventually disappear- 
ing for jjL > 1.0. We have examined such self-oscillatory 
behavior for various values of K, r, and /x, and observed 
that the oscillation frequency /o decreases with both in- 
ertia and delay, as shown in Fig. 2 for r — 1.8, 2.0, and 
2.2. 

The resulting phase boundaries in the plane of (/it, r) , 
discriminating the oscillatory state and the stationary 
one, are displayed in Fig. 3, where the coexistence of 
the two states is also identified in the region between 
solid squares and open ones. In the coexistence region, 
there exist a multitude of solutions, giving stationary or 
oscillatory behavior depending on the initial conditions. 
To explore this, we examine the behavior of the order 
parameter A = A^^^| exp(i(/)j)| with the inertia. Fig- 
ure 4(a) shows the order parameter computed via numer- 
ical simulations, as the inertia is varied between 0.1 and 
1.2 with the increment of (5/i = 0.01. We have first de- 
creased the inertia from 1.2 to 0.1, which yields the step- 
like structure consisting of the solid circles in Fig. 4(a). 
Upon increasing the inertia from given value, however, 
the order parameter does not reverse the path along the 
solid circles; instead it exhibits continuous dependence 
on the inertia, eventually following the open circles. In 
this manner a multitude of solutions or "bands" have 
been found, as shown in Fig. 4(a). Note that there are 
two kinds in shape, flat bands and bent ones. As the 
behavior of the phase velocity is examined, the system 
is stationary in the flat bands, existing for small inertia 
(/i ^ 0.3) and for large inertia (/^ > 1.1). In the bent 
bands (for 0.3 ^ /i ^ 0.9), on the other hand, the sys- 
tem displays phase oscillations. Accordingly, in the range 
0.9 /i 1.1, where two kinds of bands coexist, either 
the oscillatory phase or the stationary one is expected to 
appear depending on the initial conditions. This coex- 
istence range of inertia (for r = 2) is indeed consistent 
with Fig. 3. 

To understand analytically these numerical results, we 
first divide the population of the oscillators into two 
groups: the synchronized group {S) and the desynchro- 
nized one (D), and take the ansatz = cj)^ + uiit for D 



and 



Asm{flt + a) 



(3) 



for S. In view of the symmetry of g{uj), we consider the 
dc component 0^ to be constant. Results of numerical 
simulations also manifest that near the phase boundary 
the ac component is much smaller than the dc compo- 
nent, allowing the expansion in terms of the amplitude 
A. Upon substituting Eq. (||) into Eq. (|l|) , we thus obtain 
for the oscillators in S, to the order of A^, 
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A sin(f7t-l-a) 



^+nT. cos(0°-(/.°)sinf]T 



Acos(nt+a), (4) 



where the summations run over S. Oscillators in D, 
the fraction of which over the whole is negligibly small, 
Q^Q-{KA) /2a -j^ i^Yn out not to contribute to the har- 
monics in Q. A similar expansion for D leads to trivial 
equations. 

The zeroth harmonics in Eq. (0) leads to 



1 

j 



sm 
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or, in terms of the order parameter, 

2A^ sm^{nT/2)] 



= 1 - ^ 1 
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where is the variance of the distribution g{uj) and the 
thermodynamic limit {N — > oo) together with the strong 
coupling limit {K 3> Kc) has been taken. It is shown in 
Eq. (^) that the appearance of oscillation (A ^ 0) reduces 
A, suppressing synchronization. 

We next consider the first harmonics in Eq. (^), and 
obtain H 



XA2(1- cosier) 
isTA^ sin f7r 
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which yields only the trivial solution 57 = in the ab- 
sence of either inertia or delay. It is thus concluded that 
the (sustained) spontaneous oscillation requires both in- 
ertia and delay. Elimination of A in Eq. (Q) leads to the 
equation for the oscillating frequency f2: 
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il = tan — , 
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the nonzero solution (il ^ 0) of which may be searched 
in the range tt < fir < 27r. Equation (|^) together with 
Eq. (^) gives A^ in the form 
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Note in Eq. (||) that A approaches unity as a/K is re- 
duced to zero {(t/K — > 0). In particular Eq. gives 
the range < A^ < 1 - ctV^^, where the upper bound 
of A^ determines the boundary between stationary and 
oscillatory states. Accordingly, on the boundary Eq. (H) 
reads 
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which in turn yields the boundary 
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TT — arctan ^ / 2Kfi ^1 — 



solid squares in Fig. 3 have no counterpart in the analyt- 
ical results shown in Fig. 6. 

In summary, we have studied the synchronization phe- 
nomena in a system of coupled oscillators, each possess- 
ing finite inertia and interacting with others via time- 
delayed interactions. The interplay between inertia and 
time delay has been investigated in the temporal behav- 
ior of the phase velocity, which reveals the emergence of 
(sustained) spontaneous oscillation in the absence of ex- 
ternal periodic driving. We have also obtained the phase 
diagram, which displays oscillatory and stationary states 
in the appropriate regions of the three-dimensional space 
consisting of the inertia, delay, and coupling strength. 

We acknowledge the hospitality from Korea Institute 
for Advanced Study, where part of this work was per- 
formed, and the partial support from the BK21 Project. 
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It can be observed that the oscillatory state appears only 
for ^ > /ic = [21^(1 — <7^/K^)] . The critical value 
decreases, allowing larger regions of the oscillatory state, 
as the coupling strength K is raised and as the variance 
is reduced. The order parameter A obtained from 
Eq. (^ for X = 3 and r = 2 is plotted by the broken 
curve in Fig. 4 (b) , manifesting good agreement with the 
numerical result. The tiny discrepancy of O(10~^) be- 
tween the two curves can be attributed to the truncation 
of higher-order terms in Eq. (^). 

Solving Eq. (^ in the appropriate regime, we obtain 
the oscillation frequency as a function of inertia and de- 
lay. Figure 5(a) shows the behavior of the fundamental 
frequency il/27r on the plane of inertia and delay. For 
comparison, the dependence on inertia for given amount 
of delay t = 1.8, 2.0, and 2.2 is shown in Fig. 5(b), to- 
gether with the corresponding numerical data (previously 
plotted in Fig. 2). Note the good agreement between the 
analytical results, obtained from Eq. (^), and those from 
numerical simulations. It is manifested that the oscilla- 
tion frequency in general decreases with the amount of 
delay and inertia. 

The phase diagram in the three-dimensional space of 
{fi,T,K) is displayed in Fig. 6(a), where the boundary 
surface separates the oscillatory state from the station- 
ary one. Figure 6(b) exhibits the phase boundaries on 
the plane of inertia and delay for several values of the 
coupling strength. Also shown are the data obtained 
from numerical simulations for K = 3 (see Fig. 3), again 
demonstrating good agreement of the analytical results 
based on Eq. (11) with the simulation results. It is fur- 
ther observed that the region of the oscillatory state 
grows with the coupling strength. Note, however, that 
the analytical approach, focusing on the existence of a 
nonzero solution for il, is not able to discern the coex- 
istence of oscillatory and stationary phases from the os- 
cillatory phase. Thus the numerical data represented by 
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FIG. 3. Phase boundaries between the oscillatory phase 
(O) and the stationary one (S) are exhibited in the plane of 
(/i, r) for K = 3. Also displayed is the boundary separat- 
ing the coexistence region (C). The data represented by open 
and solid squares have been obtained from numerical simula- 
tions with five different configurations; error bars have been 
estimated by the range of the obtained values. 
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FIG. 1. Power spectrum of the phase velocity (in arbitrary 
units) for coupling strength K = 3, time delay r = 2, and in- 
ertia fj, = 0.4 (squares), 0.6 (circles), 0.8 (triangles). Only the 
fundamental pealcs axe displayed. Error baxs have been esti- 
mated by the standard deviation and lines axe merely guides 
to the eye. 
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FIG. 2. Oscillation frequency /o for /sT = 3 is shown to 

decrease with inertia /i at various amounts of delay: r = 1.8 
(squares), 2.0 (circles), and 2.2 (triangles). The sizes of the 
error bars estimated by the standard deviation are about the 
same as those of the symbols. 



2.3 r 



FIG. 4. (a) Behavior of the order parameter with inertia 
for if = 3 and t = 2, obtained numerically, (b) Comparison 
with the analytical results (broken curve). 
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FIG. 5. Behavior of the fundamental frequency f2/27r for 
K = 3 (a) in the (/x, r) plane; (b) versus /x for r = 1.8, 
2.0, and 2.2, obtained analytically. The data from numerical 
simulations, shown in Fig. 2, are also plotted for comparison. 





(b) 

FIG. 6. Phase diagram (a) in the three-dimensional 
(/Lt, T, K) space; (b) in the (/i, r) plane for several values of 
the coupling strength, obtained analytically. The data from 
numerical simulations for if = 3, shown in Fig. 3, are also 
plotted for comparison. 
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